Protocol for universal gates in optimally biased superconducting qubits 
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We present a new experimental protocol for performing universal gates in a register of supercon- 
ducting qubits coupled by fixed on-chip linear reactances. The qubits have fixed, detuned Larmor 
frequencies and can remain, during the entire gate operation, biased at their optimal working point 
where decoherence due to fluctuations in control parameters is suppressed to first order. Two-qubit 
gates are performed by simultaneously irradiating two qubits at their respective Larmor frequencies 
with appropriate amplitude and phase, while one-qubit gates are performed by the usual single-qubit 
irradiation pulses. 



Single quantum bits displaying coherence in the time 
domain have now been implemented in various super- 
conducting integrated electrical circuits Microwave 
spectroscopy 0, coherent temporal oscillations Q, and a 
conditional gate operation^] have been reported in ex- 
periments on pairs of capacitively coupled qubits. In all 
these implementations, decoherence is by far the largest 
obstacle to be overcome for applications to quantum in- 
formation processing. Yet, as Vion et al. have demon- 
strated, appropriate symmetries in circuit architecture 
and bias conditions can be exploited for suppressing to 
first order decoherence due to fluctuations in control pa- 
rameters. 

The schemes for performing two-qubit gates proposed 
so far require dynamical tuning of either the qubit tran- 
sition frequencies or of the subcircuit controlling the 
qubit-qubit interactional . The former typically requires 
dc pulses that move the qubits away from their opti- 
mal bias points for coherence, while the latter requires 
additional control lines and non-linear elements that in- 
evitably introduce additional couplings to uncontrolled 
degrees of freedom in the environment. In this Letter, 
we present a novel scheme which minimizes decoherence 
by maintaining both qubits at their optimal bias points, 
and by employing only noise-free fixed linear coupling 
reactances. Furthermore, this scheme takes advantage 
of the spread in circuit parameters occuring naturally in 
fabrication, instead of suffering from it. 

Our strategy consists of constructing a circuit with 
fixed, detuned Larmor frequencies and fixed coupling 
strengths — a sort of "artificial molecule" — and realizing 
gates with protocols inspired by those of nuclear mag- 
netic resonance(NMR) quantum computationQ. The es- 
sential difference between our "molecules" and those used 
in NMR resides in the form of the qubit-qubit couplings 
and the way they are exploited. In NMR, the secular 
terms in the coupling Hamiltonian (i.e. those that com- 
mute with the Zeeman Hamiltonian and thus act to first 
order) dominate the spin-spin interaction. Two-qubit 
gates are realized as the spins precess freely, while refo- 
cusing pulses are applied in order to "do nothing" . In our 
scheme, the coupling is purely non-secular, and has no 
effect to first order. So unlike in NMR, we must construct 
pulses to enhance the second-order effect of the coupling. 
We refer to this strategy with the (NMR-style!) nick- 



name: FLICFORQ, for Fixed Linear Couplings between 
Fixed Off-Resonant Qubits. 

The superconducting register we have in mind may 
consist of charge qubits (controlled via charges on gate 
capacitors) interacting through on-chip capacitors or of 
flux qubits (controlled via fluxes through superconduct- 
ing loops) interacting through mutual inductances. We 
focus for the moment on two-qubit registers (Fig. 1), the 
simplest that allow the realization of a universal set of 
quantum gates, leaving the extension to larger systems 
to the discussion. The optimal bias conditions for the 
circuits shown are Nf = AT| = 1/2 for charge qubits, 
where N 9 — C 9 U/2e is the dimensionless gate charge, or 
Nf = N$ = 1/2 for flux qubits, where N<t> = <P ext /<P is 
the flux frustration. Under these conditions, the systems 
become immune, to first order, to variations in the con- 
trol parameters, such as 1// charge noise in the Joseph- 
son junctions or substrate or noise due to the motion of 
trapped fluxQ. 



rH-L-Eh 

-iyi- 

- | readout 2} 




r| write 2 ^ 




-j write 1 h 






-kuhj- 



] readout TV ' — | readout 2"] — 




— | readout T] — 



FIG. 1: Superconducting two-qubit circuits for performing 
universal quantum gates at optimal bias point, (a) Charge 
qubits (Saclay style) coupled by a capacitor, (b) Flux qubits 
(Delft style) coupled by a mutual inductance. 

At optimal bias, these two-qubit systems are described 
by the reduced Hamiltonian 

H/fi = - [wf of + 2[uf (t) cos(o;[ / t) + uf(t) sin^^K 



+wfcr| + 2[ojf(t) cos^t) + oj y 2 (t) wa.{w r 2 f t)]o% 

+ , ,XX —.X —X 1 
u cr 1 a 2 \, 

where u>l/2n (wf/27r) is the Larmor frequency of qubit 



2 



1(2); uf-y I'l'K (0*2 is the frequency of the signal ap- 

plied to the "write" port of qubit 1(2); ojf (oj 2 ) and 
lo\ (w|) are the amplitudes of the in-phase and quadra- 
ture components of the applied signals, respectively, and, 
when divided by 27r, are directly interpretable as Rabi 
frequencies; and ui xx /2ir = (^t swa P^- 1 [ s the "swap" fre- 
quency (if only the crfa 2 term were present in Ti., the 
time needed to go from a product state to a maximally 
entangled state would be t swap /4). The Larmor frequen- 
cies are detuned from one another, as occurs naturally 
during fabrication, by 5 — u>f — wf, and remain fixed 
throughout. The swap frequency is fixed at the time of 
circuit fabrication and should satisfy uj xx <C S to avoid 
significant entanglement of the qubits in the absence of 
external signals. For concreteness, we consider the case 
u, xx = 0.15 = 0.01cj o , where uj = (wf + w|)/2, but 
our results do not depend critically on these values. For 
simplicity, we limit ourselves in this paper to resonant 
RF pulses constrained to obey uj^ — ojf and U} 2 = cj| 
and we perform all possible gates by playing with only 
four external knobs k/f^w^cjf and cu 2 . The difference 
between ui^ and uj r 2 suppresses cross-talk during gate 
operations, a crucial practical advantage of FLICFORQ. 

The mechanism allowing the very weak intcrqubit cou- 
pling lo xx to produce maximally entangled two-qubit 
states is easily understood in the dressed atom picture 
of quantum optics H,0|. When the RF fields and qubits 
are uncoupled, each qubit + field system has an infi- 
nite discrete ladder of doubly-degenerate energy levels, 
labelled by the qubit state |1) or |0) and the photon 
number |7V) , and separated by cu r f (Fig. 2, outer lev- 
els). Taking the qubit-field coupling into account lifts 
the degeneracy, causing the two states in each manifold 
to be split symmetrically by the field strength io\ (Fig 2, 
inner levels). The two dressed qubits may then absorb 
and emit energy at frequencies tuf ± lo\ and tu 2 ± wf , re- 
spectively. The result of the irradiation is thus to split 
the single-mode qubit line into two sidebands at these fre- 
quencies. Choosing the RF amplitudes uif 2 = 5/2 causes 
the upper sideband of qubit 1 to overlap the lower side- 
band of qubit 2, allowing the qubits to exchange photons 
of energy H(ujf — uj\) = h{u 2 + w 2 ) through the coupling 
reactance. 

The simplest protocol for generating entangled states 
in this system is to simultaneously irradiate each qubit 
with a signal of amplitude 6/2 [rJ. If we choose tUi 9 = 
6/2 and uf 2 = 0, and initialize the state to pi n — 
1 00) (00 1, the qubits will become maximally entangled af- 
ter a pulse time 4it/uj xx = 2t swa P. If the RF is then 
switched off, the system will remain in an entangled state 
until it is measured in a local basis or it undergoes deco- 
herence or relaxation. Note that this scheme allows us to 
produce entanglement on demand without dc excursions 
from the optimal bias point of either qubit. 

The rotation realized when the system is irradiated 
in this manner, which we call "D" , is not a pure erf a x 
rotation, but is rather a product of two commuting 7r/2 



rotations: 



D = (l-ia z l( T*){l + ia x aZ)/2 

= (z 1 z 2 )- 1 / 2 (x 1 x 2 ) 1 ^ 2 , 



(1) 



where we have used the rotation operator notation in 
which Xi = iaf, X\X 2 — i<jf<r 2: etc. This rotation maps 
the computational basis states to the Bell states with a 
relative phase, e.g. D: |00) -> (|00) -M| 11)) /a/2 EH, and 
can therefore be used to generate and study maximally 
entangled states. However, it is easy to verify that D 2 = 
— 1, indicating that D is a 7r rotation, and therefore is 
not universal|l2j. 

We propose to circumvent this problem by nulling out 
the <j\ a 2 factor of D. This is done by flipping the sign of 
the RF signal amplitude on one of the two qubits mid- 
way through the pulse. With this "refocusing flip" the 
unwanted erf erf rotation taking place during the first half 
of the pulse will be fully undone during the second half. 
This technique resembles the refocusing schemes used in 
NMR0 > though here we are modifying the pulse shape 
rather than performing additional 7r rotations. 

Implementing the refocusing flip leads to the pure 7r/2 
rotation (XiX 2 ) 1/2 = (1 +iafof)/^2. This rotation, 
when augmented by one-qubit tt/2 rotations, is known 
to generate the two-qubit Clifford group C 2 01 ■ So along 
with all one-qubit unitaries, (XiX^ 1 / 2 therefore consti- 
tutes a universal set of rotations. 

We can thus turn to the construction of a protocol to 
perform (/ CNOT ) the rotation corresponding to the stan- 
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FIG. 2: Energy levels of qubit + RF photons systems with (in- 
ner levels) and without (outer levels) qubit-photon coupling. 
Outer: systems have an infinite ladder of doubly-degenerate 
levels corresponding to products of a photon number state 
(green, orange) and a qubit state (red, blue). Inner: Photon- 
qubit coupling lifts degeneracy in each manifold by Rabi fre- 
quency u) v . Transitions between adjacent manifolds (wavy 
arrows) correspond to absorption/emission of a photon from 
dressed qubit system. Transition of qubit 1 at u)f — u)\ and 
qubit 2 at wf + u)\ coincide when u)\ = = 8/2, putting 
qubits on speaking terms. 
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FIG. 3: Pulse sequence for (y CN0T using FLICFORQ. (a) Am- 
plitude of in-phase (uj x ) and quadrature (w y ) components of 
resonant RF signal on first (top) and second (bottom) qubits. 
Protocol constructs (7 CN0T from sequence of five "primitive" 
7r/2 rotations. One-qubit a x and a y pulses have amplitude 
5/8 and duration t^" = 4n/S; two-qubit pulses have ampli- 
tude 5/2. The cr z rotation occurs last, where it can be ignored 
if followed by measurement in computational basis. Gate is 
completed at time t CNOT . (b) Description of sequence in 
Heisenberg picture p^l . First and last columns (red, green) 
are connected by [/ CN0T . 



dard two-qubit logical gate CNOT. We first decompose 
c7 CN0T into a sequence of rotations that draws only on 
(X1X2) 1 / 2 and one-qubit n/2 rotations. We use the se- 
quence (time runing left to right), 



xI /2 yI /2 (x 1 x 2 ) 1 / 2 y 1 



■l/2 7 l/2 



(2) 



which, as required, performs the f/ CN0T map- 
ping in the Heisenberg picture: {af, of , of , of } — > 
{of , of of > a l a 2 j a 2 Though similar in spirit to de- 

compositions of (7 CN0T given elsewhere for other sys- 
tems, e.g. expression (2) is presented in the general 
language of Pauli rotation operators, making it applica- 
ble to any physical implementation. It can, with simple 
algebra, be adapted to systems where the core two-qubit 
gate is other than (X 1 X 2 ) 1/2 ^. 

The full [7 CNOT pulse sequence is constructed by con- 
catenating the pulses generating each of the constituent 
rotations in expression (2) (see figure 3). 

We must briefly comment on how the difference be- 
tween the Larmor frequencies is dealt with. In the ab- 
sence of irradiation the natural evolution of the system 
consists of continuous rotations of each qubit about the 



a z axis, resulting in a time-dependent phase between the 
qubits that vanishes every t"v nc = 27r/c5. This phase is 
unimportant when considering one-qubit gates, as com- 
pensatory o z rotations may be realized through simple 
waiting periods |13|. However, it must be taken into ac- 
count when doing two-qubit rotations by initiating two- 
qubit pulses only at t s ^ nc = mtl vnc for integer m, i.e. 
when the qubits are in synchrony. This condition can 
be met by using one-qubit pulse amplitudes such that 
the one-qubit a x and a v rotations last t s ^ rLC . For the 
above chosen parameters, an amplitude of <5/8 for 7r/2 
pulses is convenient (finer rotations may be generated by 
weaker pulses with the same duration) , corresponding to 
tfync _ r-Q^ asso ciated timing grid is shown with 

dashed vertical lines in figure 3. This scheme is general- 
izable to multiqubit registers (see below). 

We have simulated the pulse sequence of figure 3 by 
numerically solving a set of fifteen coupled differential 
equations describing each component of the two-qubit 
density operator [16j. The simulation technique is exact 
in the sense that it does not rely on any approximations 
or perturbative expansions of the time-dependent Hamil- 
tonian. Figure 4 shows the results of a simulation of two 
representative evolutions. The simultaneous vanishing 
of each component of the two reduced density operators 
indicates the generation of a fully entangled two-qubit 
statelH. 




FIG. 4: Evolution of sample input states during U se- 
quence. Components of the reduced density operators pi = 
Trip (row 1) and p2 = Trip (row 2). pi and p2 are plot- 
ted in reference frames rotating at u)f and wf, resp. Dashed 
vertical line denotes t CNOT . Error visible at t CNOT is due 
to Bloch-Siegert shift and effect of coupling during one-qubit 
rotations!!?!!. 



level of gate fidelity can we expect from this 
We first discuss the error in one-qubit gates. 



. This 

can be nullified altogether using dynamic de- 
schemes, as is done in NMR 13j. In the present 



What 
scheme? 

Choosing the one-qubit pulse time to be t s 2 ync means 
that the coupling term uj xx will perform a parasitic ro- 
tation in this time by an angle arccos^ 31 ^ 
rotation 
coupling ; 
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system, this would be done by performing appropriately 
timed 7r rotations about a y , which anticommutes with 
the coupling term af erf . However, for the range of prac- 
tical circuit parameters 5 > O.Icj^ and lu xx < O-Olw^, 
the one-qubit gate error rate resulting from this parasitic 
rotation is already below 10 -3 , or two orders of magni- 
tude better than the fidelity of presently available read- 
out schemes 7], and the correction is an easily dispens- 
able luxury. Also, though our simulations have used only 
square pulses, realistic pulse shapes should not cause a 
significant further loss of fidelity in the one-qubit oper- 
ations, since, as is commonplace in NMR, pulse shapes 
requiring far less bandwidth could be used[jj|]. 

The two-qubit gate error rate will likely be domi- 
nated by errors resulting from imperfect RF pulses. The 
strength of the qubit-qubit interaction depends strongly 
on the amplitude of the simultaneous RF signals, so en- 
tangling gates will be sensitive to jitter or ringing in the 
pulse amplitudes. This problem could be minimized by 
using slowly-rising pulses which require less bandwidth 
rather than trying to approximate square pulses. The 
pulses could be constructed so that intended one-qubit 
rotations are implemented during the rise time. 

Nonetheless, there is still some error present in our 
simulations of two-qubit gates, even though we have 
used ideal pulses. We have verified that this can be at- 
tributed to the counter-rotating term in the rotating wave 
frameworkp^. as the qubits are irradiated with strong 
fields for many Larmor periods during a two-qubit gate. 
This error can be reduced by choosing a stronger coupling 
uj xx , thereby reducing the time required to generate en- 



tanglement, or by reducing the detuning 8, which reduces 
the required field strength. Since these changes would 
increase the one-qubit gate error rate due to the fixed 
coupling, an NMR-style decoupling scheme will likely be 
needed once we require a gate error rate < 10~ 3 . 

We believe a main advantage of the gate scheme pre- 
sented in this paper is that it can be directly generalized 
to multiqubit registers with minimal extra hardware. A 
fixed linear coupling reactance between all pairs of qubits 
could easily be achieved by coupling each qubit to a com- 
mon superconducting island, loop or cavity. Selective 
one-qubit gates would be realized by applying an RF sig- 
nal at just the target qubit transition frequency, while 
the protocol generating D or the universal two-qubit ro- 
tation (X1X2) 1 / 2 could be realized on any pair by simul- 
taneously applying RF signals at the resonant frequencies 
of the two targeted qubits. Since each qubit in the regis- 
ter would be detuned from all the others, all these write 
pulses could be multiplexed on a single RF control line, 
a decisive advantage in seeking to limit stray couplings 
to the environment or crosstalk between qubits. Apply- 
ing D to several qubits in a pairwise fashion would al- 
low the direct production of multiqubit entangled states 
of the form \GHZ) = (|0)® n + e^|l)®")/%/2, which, 
for n > 2 can display maximal violations of Bell-type 
inequalities 0. 
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